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Abstract 

Ph 

c/j A quite straightforward approximation for the electrostatic interaction between two perfectly 

conducting surfaces suggests itself when the distance between them is much smaller than the 

q 

characteristic lengths associated to their shapes. Indeed, in the so called "proximity force ap- 

O 

proximation" the electrostatic force is evaluated by first dividing each surface into a set of small 

>, 

^ flat patches, and then adding up the forces due two opposite pairs, the contribution of which are 

Oh 

approximated as due to pairs of parallel planes. This approximation has been widely and suc- 
cessfully applied to different contexts, ranging from nuclear physics to Casimir effect calculations. 

CO 

We present here an improvement on this approximation, based on a derivative expansion for the 

o 

electrostatic energy contained between the surfaces. The results obtained could be useful to dis- 
cuss the geometric dependence of the electrostatic force, and also as a convenient benchmark for 

CN 

H numerical analyses of the tip-sample electrostatic interaction in atomic force microscopes. 

> 

i 
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I. INTRODUCTION 



A standard problem in electromagnetism is to compute the electrostatic force between 
conducting bodies, or its close relative: the calculation of the capacity of a system of con- 
ductors. The simplest example is the case of parallel plates separated by a distance much 
smaller than the characteristic size of the plates, which are held at a fixed electrostatic po- 
tential difference. Albeit not as straightforwardly as in that example, some other systems 
admit analytical exact solutions also; indeed, that is the case for two eccentric cylinders, for 
a cylinder in front of a plane, and also for a sphere in front of a plane. 

This problem is also of considerable practical relevance in Electrostatic Force Microscopy 
(EFM) and its variants p], which are based on the interaction between biased atomic force 
microscopy (AFM) tips and a sample. The same applies to the experimental determination of 
Casimir or gravitational forces between conducting bodies, as residual charges or potentials 
produce unwanted forces that must be subtracted in order to determine the sought-after 
force [2]. 

Of course the electrostatic force between bodies of arbitrary shape can, in principle, be 
computed by solving numerically the Laplace equation with adequate boundary conditions. 
However, analytic or semi-analytic methods are always welcome, as ways to improve the 
understanding of the geometric dependence of the force, and also to be used as simple 
benchmarks of fully numerical computations. For instance, in the context of EFM/AFM, 
analytical models for the tips have been developed, and additional exactly solvable models 
have been found, like the case of an hyperboloid in front of a plane j3]. Generalized image- 
charge methods have also been proposed [I], in which the tip and the sample are replaced 
by a set of fictitious charges, whose intensity and positions are found numerically. 

An interesting analytical approach was introduced by B. Derjaguin in 1934, who de- 
veloped an approximate method to compute the Van der Waals force between macroscopic 
bodies assumed to be close to each other [51 E] • The approximation assumes that the surfaces 
of the bodies can be replaced by a set of flat patches and that at short distances the domi- 
nant contributions correspond to pairs of patches (one on each surface) which are closest to 
each other. Moreover, the interaction is supposed to be additive. In this way, it is possible 
to compute the force between gently curved surfaces from the knowledge of the interaction 
energy for flat surfaces, as long as the radii of curvature of the surfaces are much larger than 
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the minimum distance between bodies and when the surface normals in opposite patches 
are almost parallel. Later on, the same idea was applied to nuclear physics, under the name 
of Proximity Force Approximation (PEA) or Proximity Force Theorem, in order to compute 
the interaction between nuclei E] • It has also been widely used to compute Casimir (or 
retarded Van der Waals) forces between neutral macroscopic objects [2J. It should be clear 
that, mutatis mutandis, the Derjaguin approximation and its ulterior developments can also 
be applied to the analysis of electrostatic forces, at least between gently curved conducting 
surfaces. 

In spite of the simplicity and usefulness of the PEA, for many years there was a stumbling 
block to its progress, since methods to asses its reliability, or to compute the next to leading 
order (NTLO) correction were lacking. This situation left, as the only alternative of assessing 
the PFA reliability, its comparison with the rather few examples for which an exact result 
was available. 

In an attempt to improve that situation, we have recently shown [Sj, in the context 
of Casimir physics, that the PFA can be considered as the first term, in an expansion in 
derivatives of the surfaces shapes, of the interaction energy. In this way, it is now possible 
to improve the PFA by computing the NTLO corrections. 

In this paper we carry out this idea in the context of electrostatics, as already suggested 
in [9] . Our aim is twofold: on the one hand, to show the potential usefulness of the improved 
PFA in computing electrostatic forces. On the other hand, we believe that this attempt to 
present the improved PFA in a simpler context, may help one to gain intuition about its 
applicability in more complex scenarios. 

II. PROXIMITY FORCE APPROXIMATION 

Let us assume, for simplicity's sake, that the system consists of a gently varying surface 
in front of a plane (both perfect conductors). The plane is at z = 0, and we assume it to be 
grounded. The curved surface is described by a single function z = ip(x,y), and is assumed 
to be at a constant electrostatic potential V. We shall use the notation x» = (x,y). The 
electrostatic energy contained between surfaces is then given by: 




(1) 
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The force between the conductors can be obtained by computing the variation of U under 
a rigid displacement of one of the surfaces, while the capacity of the system is given by 
C = 2U/V 2 . 

In the PFA, this system is replaced by a set of parallel plates (see Fig.l). The intensity 
of the electric field between parallel plates separated by a distance if) has the ^-independent 
value V/if). Therefore, the electrostatic energy is approximated by 



When computing Van der Waals, nuclear or Casimir forces, the rather involved nature of 
the interaction considered can render the interpretation of the approximations involved in 
the PFA somewhat obscure. Here, on the contrary, the physical assumption is quite clear: 
the electric field between opposite patches is regarded as constant and directed along the 
line that joins both patches. 

Eq.(|2]) can be used to estimate the electrostatic energy in many interesting situations. 
To illustrate this point, let us consider a cylinder of length L and radius R in front of a 
plane, and let us denote by a the minimum distance between surfaces. The function if), that 
describes the part of the cylindrical surface which is closer to the plane, reads: 



with — x m < x < xm < R in order to cover the part of the cylinder which is closer to the 
plane. It is reasonable to assume that xm/R = 0(1) < 1, and under this assumption the 
final result will not depend on %, as long as R 3> a, where the PFA is expected to give an 
accurate value of the electrostatic energy. Inserting Eq.Q into Eq.(|2]), the integral can be 
computed in terms of elementary functions. Expanding that result for a <C R we obtain 



which is independent of %. Note that, when the cylinder is close to the plane, the electro- 
static force is proportional to a~ 3 / 2 . 

For this geometry, the exact electrostatic interaction energy is known [TO] : 




(2) 




(3) 




(4) 
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FIG. 1. In the Derjaguin approximation, the interaction between a gently curved surface and a 
plane is approximated by that of a set of parallel plates. For each pair of parallel plates, border 
effects are ignored. 

In the limit of a/R <C 1, this exact electrostatic energy indeed reduces to the PFA result 
given in Eq.Q. For later use, and to estimate the relevance of the corrections to the PFA, 
we expand the exact result including the NTLO in a/R: 



that the NTLO produces a noticeable improvement of the PFA. We will now compute 
such improvement for an arbitrary geometry, of which the exact electrostatic energy in not 
necessarily known. 

III. IMPROVING THE PFA: DERIVATIVE EXPANSION 

In order to improve the PFA, we shall use the fact that the electrostatic energy may be 
thought of as a functional of the shape of the surface. This functional will be, in general, 
a non-local object, which becomes local when the surfaces are sufficiently close and parallel 
to each other. To interpolate between those two situations, we shall therefore assume that 
the electrostatic energy can be expanded in derivatives of if). Indeed, one can think of the 
condition \Vip\ < 1, as the existence of a small, dimensionless parameter. This translates 
into the fact that the upper surface is almost parallel to the plane. 

Including terms with up to two derivatives, the electrostatic energy has to be of the form: 





In Fig.2 we plot the ratios U cp /Up P FA and U cp /U{ 



rCp 

NTLO 



function of a/R. We see 




(7) 
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FIG. 2. Ratio of the exact energy and PFA (upper curve) for the cylinder in front of a plane, as 
a function of a/R. For the same geometry, the lower curve corresponds to the ratio of the exact 
energy and the approximation that includes the NTLO correction. The inclusion of the NTLO 
improves notably the approximation. 

for some functions V e s and Z [11]. The result must be proportional to e^V 2 , and must 
reproduce £/pfa f° r constant ip. Moreover, as there are no other dimensional quantities 
beyond ip, dimensional analysis imply that both functions V c r and Z must be proportional 

to if)' 1 : 



where /3 em is the only (numerical) coefficient to be determined. Note that this coefficient is 
independent of the nature of the surface being considered; therefore, it may be obtained once 
and for all from its evaluation in a single case. Indeed, in order to obtain the coefficient /3 em , 
the simplest approach, when an exact solution of the problem is available, would be to read 
its value by expanding the exact solution. For instance, for the particular configuration of 
a cylinder in front of a plane, one can insert Eq.(|3| into Eq.(|8]), perform the integrals, and 
expand the result in powers of a/R. Comparing this with the expansion of the exact result 
given in Eq.([6]), it is straightforward to show that (3 cm = 1/3. Of course the same result 
could have been obtained from any other particular example. 

Nevertheless, the coefficient (3 em can also be computed from first principles, as described 
in the next section. 




(8) 
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IV. SOLVING THE LAPLACE EQUATION 



It is perhaps more satisfying to prove that /3 em = 1/3 without assuming a particular 
shape for the curved surface, but rather for a rather general class of surfaces. To that end it 
will be sufficient to derive it by assuming that ip{x\\) = a + i](x\\) where a is the mean value 
of ip, and <C 1. For this class of surfaces, the derivative expansion of the electrostatic 

energy reads 



e V 



DE 



2a 



d 2 x 



(9) 



1+Q) +P em \V V \ 2 + 0(r ] 3 ) 

and therefore he coefficient (3 em can be read from the quadratic dependence of Ube on the 
derivatives of rj. 

The electric field is given by E = — V0, where </> is the electrostatic potential, which 
satisfies the Laplace equation with the appropriate boundary conditions on the surfaces: 

V72, 



0. 



>z=Q 



0. 



>Z=1p 



V. 



(10) 



To proceed, we follow the approach of [12] . to trade a simple boundary condition on 
a complex surface into a complex boundary condition on a flat surface. We expand the 
boundary condition on the upper surface as 

1 



V = <j>(x\\,a + 7)) = 4>(x\\,a) +r)d z <j)(x\\,a) + -rj d z <j>(x\\,a) + ... 
and we expand the electrostatic potential as 

4> = 4>o + 4>i + 4>2 + ■ ■ ■ , 



(11) 



(12) 



where 0o is the solution for parallel plates separated by a distance a, that is <po = Vz/a, 
and (pi is assumed to be 0{rf). 



Inserting Eq.(12) into Eq.(ll), we see that (f>i and <p2 must satisfy the following boundary 



conditions on z = a: 



(j)i{x\\,a) 



V 



■q(x\\ 



(f) 2 (x\\,a) = -r](x\\)d z (t)i(x h a) 



(13) 



Thus, the initial problem has been converted into a problem of parallel plates for <pi and 02, 
with more involved boundary conditions. Up to quadratic order, the electrostatic energy 
reads 



DE 



d 2 x 



a+r] 



dz 



V 2 V V 

— + |V0!| 2 + 2-0^ + 2-0^2 



(14) 



It is easy to check that the last two terms in Eq.(14) do cancel: 

ra+T] 



d A x\\ I dz (d z (f)i + d z <p2) 
o 



d x\\ r)(x 



(15) 



where we have used that 



+ 



>2)\a+ V 



rj + o(rj 2 ) and that the mean value of rj vanishes. 



Therefore the correction to the PEA is given solely by the term proportional to |V0i 



The order-77 solution <f>i is given by 



h{x\\,z) = -- 



V 



d 2 h ik„r„~r, .sinhkz 



(16) 



a J (27r) 2 /v " y sinh /ca 
Indeed, it is simple to check that it satisfies both the Laplace equation and the correct 
boundary condition (we are using the notation k = \k»\). The solution for 2 is similar- 
looking, although as already pointed out its explicit expression will not be necessary in 
order to evaluate the second-order energy. 



Using the solution Eq. ( 16 ) we obtain: 



d*x\V(j)i\ 2 



a 2 



d 2 k\\ 



-m\\)\ 2 nh 



where 



a , e 

dz 2 — I 

sinh ka 



1 + 2 sinh 2 kz) 



1 



ak 2 
~3~ 



(17) 



(18) 



(note that the integration in z runs up to a instead of a + rj since the integrand is already 
of order rf). Therefore 

e V 2 



d 6 x\V(f)j 



d x 



2 J 1 r 1 2a 
which is of the form of Eq.([9]), and shows that /3 em = 1/3. The final result is then 



e V 2 



DE 



d 2 x\\ 



(19) 



(20) 



as expected from the example of the cylinder in front of a plane. Eq.(20) is the main result 



of this paper. The vertical component of the force on the upper surface can be computed 
by taking the variation of the energy under an infinitesimal translation ip — >■ ip + 8, that is 



dU 



DE 



e V 2 



d x 



II ^2 



(21) 



dS 2 J " ~" i) 2 

Note that, as the conductors are not isolated (they are kept at a fixed potential difference), 
the force is given by the derivative of the energy with respect to the position, with a positive 
sign. 
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V. TIP-SAMPLE ELECTROSTATIC INTERACTION 



In this section we use the improved PFA to analyze the problem of the tip-sample elec- 
trostatic force in an AFM. Typical geometries to model the interaction are sphere-plane, 
paraboloid-plane, hyperboloid-plane, and sphere-ended cone in front of a plane, among oth- 
ers [3]. In all these cases, it is possible to obtain the electrostatic force using the improved 
PFA. The results will be accurate when the tip is sufficiently close to the sample. 



A. Sphere-plane 

As a first example, we describe the case of a sphere of radius R, in front of a plane at 
a minimum distance a (a <C R)- The surface of the sphere cannot be described by a single 
valued function z = ip{x\\). We will just consider the region of the sphere which is closer to 
the plane, and we will see that it still give quantitatively adequate results even beyond the 
lowest order. The function ip is 



^p) = a + R[l- X ll-^}, (22) 

where we are using polar coordinates (p,4>) in the plane {x,y). This function describes the 
hemisphere when < p < R. The derivative expansion will be well defined if we restrict 
the integrations to the region < p < pu < R- As in the example of the cylinder in 
front of a plane, described in Section II, we will assume that Pm/R = 0(1) < 1. Inserting 



this expression for ip into the derivative expansion for the electrostatic energy Eq.(20), one 
can perform explicitly the integrations and obtain an analytic expression Ube{pm, R)- 
Expanding the result for a/R<^i 1 we obtain 

u DE ~ U? FA + /i + ~ (^ P FA + h) , (23) 

where the leading contribution is given by 

Up FA = -nV 2 e R \og{2a/R) , (24) 

and f\ t 2 are functions of Pm/R. Being independent of a, fx can be discarded if one is 
interested in computing the force between the sphere and the plane. The term proportional 
to f2 gives a next to NTLO contribution to the force. Therefore we write the final result as 

%E-^ FA (l + ^). (25) 



Note that, as long as a <C R, the force will not depend on pu'- 

cIUde nV 2 e R ( la /a 



1 + • (26) 



da a \ 3R \R; 

This result could also be obtained directly from Eq.(21), avoiding the appearance of the 
PAf-dependent constant terms f 1>2 . 

As already mentioned, this geometric configuration admits an exact solution [14]; at short 



distances that exact result reduces to Eqs.(25) and (26), as discussed in Ref. [15j (we include 
some details in the Appendix). 

Finally, it is worth to remark that at short distances the electrostatic force for this 
configuration is proportional to a -1 , milder than the cylinder-plane geometry (a~ 3//2 ) and 
than the parallel plates (a~ 2 ). This is of course a consequence of the different effective areas 
of interaction in each case, and has its counterpart in Casimir physics [13J. It is a simple 
property that could be emphasized in elementary discussions of electrostatics. 



B. Paraboloid-plane 

For the case of a paraboloid-like tip, the surface is described by 

^(P)=a+^ (27) 

with < p < pm < R- k — R~ l is the local curvature at the bottom of the paraboloid. The 
derivative expansion for the electrostatic interaction energy is 

^^^5a(i-~). (28) 

It is not surprising that the leading order coincides with that of the sphere-plane interaction, 
since both surfaces have the same curvature at the point closest to the plane. However, the 
NTLO corrections are different, and in fact have opposite signs. These corrections can 
therefore be used as benchmarks of numerical calculations, since the short distance regime 
is typically the complex from a computational point of view. 

C. Hyperboloid-plane 

In this case the surface is described by the function 



^(p)=a-R 1 + R 1 Jl + ^, (29) 
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where the curvature k at the bottom of the hyperboloid is k = R^ 1 , with R = R\/R\, and 
R1/R2 is its asymptotic slope. 

As in the previous examples, the derivative expansion of the electrostatic energy can be 
computed in terms of elementary functions. The result is 

,2 i? 2 cT 



U B e ^ ^pfa 



l 3 R\'R\ 



(30) 



It is remarkable that sign of the NTLO correction depends on the parameters of the hyper- 



boloid, and vanishes for the particular asymptotic slope R1/R2 = a/3/2. 

VI. GENERALIZATIONS 

The results presented so far can be generalized in several directions: 

A. Beyond the next to leading order 

It is in principle possible to extend the results in order to include more derivatives of the 
shape function. For example, coming back to the case of a surface in front of a plane, the 
improved PFA would be, including the next to NTLO corrections: 

U DE = e -^fj rf 2 xn^(l + i|V^| 2 + 7i|V^| 4 + 7 2 ^|V^| 2 V 2 ^ 

+ 73^ 2 V 2 ^V 2 ^ + ^ 2 d Q d^d a d^ + 7 5?A 3 V 2 V 2 V0, (31) 

where 71, 72, ...75, are numerical coefficients that can be determined extending the calculation 
of Section IV up to 0{rj 4 ) 



The structure of the improved PFA is useful to analyze its range of validity: the NTLO 
corrections are small when |VV>| <C 1 or, in other words, when the curved surface is almost 
parallel to the plane. Higher order corrections will be negligible when, in addition to this 
condition, the scale of variation of the shape of the surface is much larger than the local 
distance between surfaces. 



B. Flat surface coated with a dielectric layer 



The improved PFA can be extended to the case in which the sample (flat surface) is a 
perfect conductor coated with a layer of width d of a material characterized with a constant 
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permittivity e, as considered in the context of EFM [T7]. Dimensional analysis implies that 
the derivative expansion for the electrostatic energy will be of the form 

t/ DE * ^ / rf 2 x,| - -) + /3(-, -)|V^| 2 1 , (32) 

2 J if) I e if) e if) 

for some functions a and /3. We do not dwell here with the determination of the function 

(3, but just want to remark that a can be obtained by elementary considerations for a 

parallel-plate capacitor with a dielectric layer: 

. e d. 1 , . 

a(- -) = 33 

C. Two gently curved surfaces 

For the case of the Casimir energy, the result of Ref. [9] has been generalized to the case 
of two gently curved surfaces in Ref. [18J. This generalization can also be applied to the 
electrostatic case. Let us assume that the surfaces are described by two functions ipi(x\\) 
and ^2(^11) that measure the height of the surfaces with respect to a reference plane E. The 
derivative expansion for the electrostatic energy reads, in this case: 

£^1,-02] = ^ / ^ x \\^ [1 + /3i|V^i| 2 

+ /3 2 |V^ 2 r + Px W>i ■ W> 2 + P- z • VV>i x V^ 2 ) + ■■■], (34) 

where if> — \ip2 — ^i| is the height difference, and dots denote higher derivative terms. 

As the energy must be invariant under the interchange of ip\ and 1P2, Pi = P2 and /3_ = 0. 
Moreover, in order to reproduce the previous result we must have /3i = /32 = 1/3. Finally, 
the coefficient [3 X can be determined using the fact that the energy should be invariant under 
a tilt of the reference surface £ [IB] or, alternatively, under a simultaneous rotation of both 
surfaces. For an infinitesimal rotation of each surface with an angle e in the plane (x, z) we 
have 5ipi = e(x + if>id x if)i) , for i = 1,2. In order to simplify the calculation we can assume 
that, initially, ip\ = and that if>2 is only a function of x. Computing explicity the variation 
of Ude to linear order in e one can show that 

SU DE = (3 X = \ , (35) 



and therefore 

£ / de[V , i,'02] = ^-^ir- I d' 2 x\ 



^V 2 f l2 1 



3 



l + ^(|V^i| 2 + |V^2| 2 + V^-V^ 2 ) 
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(36) 



While the equality 0\ = 02 is valid for any interaction (as long as the surfaces are 
identical), the fact that (3 X = Pi is valid only for the particular case of the electrostatic 
interaction, in which the leading term is proportional to ip~ l (i.e. is not valid for the 
Casimir energy, that is proportional to for parallel plates). 



Computing the variation of the electrostatic energy Eq.(34) with respect to translations 
or rotations of one of the surfaces, it is possible to obtain the vertical and lateral components 
of the force, as well as the torque produced by the presence of the other surface. 



VII. FINAL REMARKS 



The PFA is a very useful approximation introduced by Derjaguin many years ago to 
compute Van der Waals forces between macroscopic bodies. The approximation has been 
subsequently widely used in rather different contexts, like nuclear physics and the Casimir 
effect. Very recently, the approximation has been interpreted as the leading order in a 
derivative expansion of the energy, and this interpretation allowed to compute, for the first 
time, the NTLO corrections in a systematic way. In this paper we have pointed out that the 
improved PFA is a powerful tool which may also be used to compute the electrostatic inter- 
action between conducting surfaces. Indeed, with this improved approximation, it is possible 
to compute the electrostatic force with high accuracy when the surfaces are sufficiently close 
to each other. 

Moreover, given the simplicity of the expressions for the energy and the force, the NTLO 
could be used to test numerical methods aimed at a calculation of the force for arbitrary 
surfaces. 



NOTE ADDED: 

While this paper was under review, the revised version of Ref.[T8] was published |19j . 
In the new version, the authors mentioned that they have successfully applied the gradient 
expansion to analyze the electrostatic interaction between a sphere and a plane, although 
there are no details of the calculations. 
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APPENDIX 

The electrostatic problem of a sphere of radius R in front of a plane can be solved exactly: 
if the potential difference between the plane and the sphere is V, the electrostatic energy is 
given by jH] 

U sp = 2ire RV 2 sinh T V , * _. , (37) 

t-^ sinh (nl 

n>l v ' 

where coshT = 1 + a/R. In order to obtain an analytic expression in the limit r — > we 
write 

S = Y - l — = S - [°° ^— +^ln(cothr) 
smh(nr) J l smh(nr) V 

= ^ + ^ln(cothr) + 0(r), (38) 



where 7 = 0.5772. Replacing this expression into Eq.(37) and expanding the result for small 
a/R we obtain 



l/sp ~ -ne aV 2 \n(2a/R) jl + ~ + ^ 



a/R 



(39) 



Jn(a/i?), 

where we omited an irrelevant constant term. This result coincides, up to the NTLO, with 
the one obtained using the derivative expansion Eq.(25). 
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